Sub-Dominant Pairing Channels in Unconventional Superconductors: 

Ginzburg-Landau Theory 
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A Ginzburg-Landau theory is developed for unconventional superconductors with the three relevant 
singlet pairing channels (i.e., s-, d^2_y2- and dxy channels). Various consequences of the sub- 
dominant channels (i.e., s- and da^^-channels) are examined in detail. (1) In the case of a d^2_y2 +is- 
wave superconductor, we reproduce an earlier result that there is a second order zero-field transition 
from the pure d^2_y2-phase to the time-reversal-symmetry- (T-) breaking d^2_y2 -I- is-phase at the 
temperature Tds- The structure of a single vortex above and below Tds is four-fold and two-fold 
symmetric, respectively. (2) In the case of a d^2_y2 +idxy-wa.ve superconductor, there is also a second 
order zero-field phase transition from the pure d^2_y2-phase to the T-breaking d^2_y2 + idxy-wa.ve 
phase at the temperature Tr)r>i . In contrast to the case in a d^2_y2 + is-wave superconductor, the 
subdominant phase can not be induced by vortices above T^u^i. Below the T-breaking transition, 
the sub-dominant phase in the mixed state is nontrivial: it survives at low fields, but may disappear 
above a field (increasing with decreasing temperature) presumably via a first-order transition. (3)By 
including the strong coupling effects, a T-breaking coupling term between the d^2_y2- and d^^^-waves 
is found to have significant effects on the low temperature behavior of d^2_y2 +idxy superconductors. 
In a magnetic field, a d^2_y2 -\- id^y state is always established, but the field-dependence of d^y- 
amplitude above Tjjjji is different from that below Tjjjji. Above but not very close to Tjjjji, the 
induced minimum gap Aq oc B/{T — Tub')- 
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I. INTRODUCTION 

The phase-sensitive experiment of the unusual Joseph- 
son effect in high-Tc superconductors evidenced thft 
d3;2_j^2-wave symmetry of Cooper pair wave functions.cl 
Thus at least the problem as what is the dominant 
channel for electrons to pair in high-Tc superconduc- 
tors has been resolved. However, a new puzzle arises as 
what is the sub-dominant pairing channel, which seems 
to have manifested in a number of experiments, e.g., 
the observation of surface-induced broken time-reversal- 
symmetry (T hereafter) in YBCO tunnel junctions,a the 
observation of fractional vortices trapped in a boundary 
j unction J3 and the abnormal field dependence of thaJow 
temperature thermal conductivity Kg in BSCCO,qEI to 
name only a few. At surfaces, the sub-dominant chan- 
nel can show up because of a suppression of the effec- 
tive critical temperature in the dominant d2;2_y2-channeL 
Rather differently, the anomalous Ke reported in Ref.cl 
might Keciuire a mechanism to produce hulk T-breaking 
states. EIu This scena|do might work above bK in the sam- 
ples reported in Ref.ffl but failed to explain the data at 
sub-Kelvin temperatures in Ref.El. Instead, field-induced 
quasiparticle excitations with pure (ia;B]_y2-wave gap func- 
tion explain qualitatively these dataJs 

As a modeling study, we consider the relevant sin- 
glet sub-dominant pairing channels in an unconventional 
d- superconductor to be the s- and d^j^-channels, which 



may be relevant to high temperature superconductors. 
For attractive pairing interactions in the s- and dx2-y2- 
channels, or in the dx2-y2- and (ia;j,-channels, the the 
uniform superconducting state in the absence of mag- 
netic field has been shown to be either a pure state in 
the dominant channel, or a state with a broken ivpe- 
reversal-symmetry involving the relevant channels. Op In 
order to study the vortex states, we develop a Ginzburg- 
Landau (GL) theory that includes the above-mentioned 
sub-dominant channels. We discuss the consequences of 
these channels on the properties of the superconductors, 
such as the field dependence of the thermal conductivity. 
The structure of this paper is as follows. The microscopic 
derivation of the GL theory (in the weak coupling limit) 
is presented in section II. The properties of dj;2_y2 -\-is su- 
perconductors are described in section III. The properties 
oid^2_y2 -\-idxy superconductors are addressed in section 
IV. A phenomenological strong coupling correction to the 
weak coupling theory is developed and discussed in sec- 
tion V with respect to the d^2_y2 -\-idxy superconductors. 
Section VI contains a summary and some concluding re- 
marks. 



II. DERIVATION OF THE GINZBURG-LANDAU 
THEORY IN THE WEAK COUPLING LIMIT 

The simplest way to derive microscopically a 
Gingzburg-Landau theory for superconductors is to con- 
sider the Bardeen-Cooper-Shrieffer (BCS) gap equation. 
At zero center-of-mass momentum for the Cooper pairs, 
the gap function is 
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which follows routinely from the effective BCS Hamil- 
tonian. Here Ef^ = yjej. + jA^p is the quasi-particle 
excitation spectrum, and /3 = 1/T is the inverse tem- 
perature. (We are using the units fc^ = ?i = c = 1, 
and are measuring energies from the Fermi surface.) In 
the weak coupling limit the summation over momentum 
in the gap equation is understood to be restricted near 
the Fermi surface. Before going into further details, we 
need to specify the pairing interaction Vk^k' ■ Due to the 
tetragonal symmetry of the copper planes, it is rather 
reasonable to assume the following form for Vk_k' : 

VkM' ^Vs+ V,.^y2{kl - kl){k'l - k'l) 

' ^* xy'^x'^y'^ x^ y 



= Vs + Vx^-y^ cos 29k COS 26k' 
+ Vxy sin 26k sin 29k' ■ 



(2) 



Henceforth the subscripts s, x^-y^ and xy are channel in- 
dices referring to the s-, dx2_y2- and d^^-channels of the 
pairing interaction, respectively. 6k is the angle between 
the wave vector k and the a-axis in the copper-oxide 
plane. By symmetry, the gap function can be expressed 



l^k = S + D cos 29k + D' sin 29k , 



(3) 



where S, D and D' are complex amplitudes of the gap 
function in the s-, dx^-y^- and d^^^- channels, respectively, 
and are also called as the order parameters. Inserting 
Eqs. (H) and (||) into Eq. (|l]), expanding the gap equa- 
tion to the third order in the order parameters [using the 
Poisson identity {l/2E)ta.nh{PE/2) = TY,JloI + E'^)~^ 
with LJn being the Fermion Matsubara frequency] , we ob- 
tain the GL equations in the uniform case, 

= -N{Q)Vs {2asS + 7 [2\S\'^S + \D\'^S + \D'\'^S 
+ {S* D + c.c)D + {S* D' + c.c.)D']} ■ 

3, 
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N{Q)Vxy\aD'D' + ^ 



hD'\^D' + -\D\^D' 



\S\^D' + {S*D' + c.cjS + -{D*D' + c.c.)D 



Here N{Q) is the normal state density of states (DOS) 
at the Fermi surface, Ui — InT/Ti {i — S,D,D') 
with Ti the bare critical temperatures at the i- 
channel: Ts w 1.14wcexp(-l/7V(0)K), and Td^d' ~ 
1.14uJcexp{—2/N{0)Vx2_y2j.y) ( with lUc the energy 
cut-off for the pairing interaction). Finally 7 = 
/^_2 deTJ2ni^n + ^'^)~'^- On the other hand, the rigidity 
of the order parameters can be probed by the long wave- 
length (^ Xp) spatial variation of the order parameters. 
For this purpose, it suffices to consider the linearized gap 
equation but with a small center-of-mass momentum q 
(<C kp) for each Cooper pair. 
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where k± — k ± q/2, and /(e) is the Fermi distribu- 



tion function. We now assume A 



k,q 



Sq + D„ COS 26k 



D' sin 29k , and expand the above equation to the second 
order in q to find that the Fourier components Sq, Dq 
and D' satisfy 



<ll)Dq 



= -N{Q)V, {2asSq + K [2q^Sq + {ql 

+ iqxqy + qyqx)Dq]} ; 

= -N{0)Vx2^y2 [aoDq + K [q^Dq + {ql - ql)Sq\ } ; 
= -'N{Q)Vxy {ao'D'g + K [q^D'g + {qxqy + %fe)^,] } , 

where K = Vpj/8 refiects the rigidity of the order pa- 
rameters with vp being the Fermi velocity. Resorting to 
real space, and using the gauge invariant gradient oper- 
ator n = — iV — 2eA (with A the vector potential) in 
substitution of — «V, and finally combining the nonlin- 
ear homogeneous tems in Eq. (Q) , we obtain the complete 
GL equations (up to third order in the order parame- 
ters) describing pairing in all of the three most relevant 
spin-singlet channels. 
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+ {S*D + c.c.)S + j{D*D' + c.c.)D' 


+K [tfD + (H^ - Wl)S] = 0; 
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h(S*D' + c.c.)S + {D*D' + c.c.)D] 
-K [li^D' + (H.Hj, + nyIVx)S] = 0. 



For completeness and general purposes, we need a GL 
freie energy functional. Following from Eqs. (|g), it is given 
byl 
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where uniformity along the c-axis is assumed, and /^^ 
denotes integration over the a6-plane. The prefactor of 
-/V(0)/2 could be obtained from a microscopic derivation 
of the super-current, but it can also be obtained by the 
inspection that the free energy reduces to the supercon- 
ducting ground state energy at zero temperature. The 
last term in F is the magnetic energy. The equilibrium 
super-current can now be easily obtained from the func- 
tional derivative of the free energy with respect to the 
vector potential, SF/SA = 0. It can be seen clearly that 
the GL equations (or the free energy functional) would 
be identical to that derived previously for dx2-y2 + is 
superconductors if the d^^y-channel were shut-down, and 
would be equivalent to that investigated for dx2_y2. +idxy 
superconductors if the s-wave channel were disregarded. 
Moreover, the gradient terms for D and D' are symmetric 
with respect to each other under a rotation of 7r/4. 



III. Dx2-y2 +/S SUPERCONDUCTORS 



This case has been studied previously,E3T3Ej so that 
we shall first recapitulate some essential points that in- 
terest us, and add some novel discussions. 

We assume that Tn > Ts as is the case in high- Tc 
superconductors. Define asD = otsl\p^D\ with ao < 
0. Optimizing the free energy Eq. (m one finds that 
the uniform bulk phase below the highest bare critical 
temperature Tjj is a pure D phase at asD > —2/3, or 
above a critical temperature Tsd — Tg/T^, at which the 
system undergoes a second-order phase transition to a 
T-broken phase (5, D) ex (±177,!), with 77 = \S\/\D\ = 

^-{3asD + 2)/{4 + 2asD)- 



Although there is no uniform 5* at T > Tsd, S can be 
induced by inhomogeneities, such as surfaces (or twin- 
boundaries), disorders and vortices, due to the mixed 
gradient terms in Eq. (|6|). Of special interest is the S*- 
component generated by vortices. At T > Tsd, the 
S'-component is sub-dominant and vanishes where D 
is uniform (e.g., far from the vortex core). Thus it 
can be obtained perturbatively. Assume D ^ e'^ in 
cylindrical coordinates and in the gauge A — A{r)d. 
Since 11^ — Ul transforms as e*^'^ and scales as l/r^ 
at r ^ 00, wc immediately see that to the first order 
S ~ a(r)e'^'^ + 5(r)e^*^ (the concrete forms of a{r) and 
b(r) are unimportant, except that both a{r) and b(r) van- 
ishes as r and 1/r^ as r ^ and 00, respectively) so that 
IS'p ^ a{r)^ + bir)'^ + 2a{r)b{r) cos 46 . (For complex a 
and 6, there are ^-independent phase shifts in the argu- 
ment of the cosine function, but they do not alter our 
general conclusion. This is also the case in similar dis- 
cussions hereafter). Therefore \S\'^ is four- fold symmet- 
ric near the vortex core and vanishes far from the core. 
This rpault.. has been obtained previously by many au- 
thors, ll3"E3o in different contjexls, and has been verified 
by our numerical simulations .EjHj In particular, Franz et 
al were able to obtain explicit analytical as rWjell as nu- 
merical solutions for the single vortex state|13 which is 
useful for a quantitative comparison between the theory 
and the experiments. The relative phase of these two 
components varies continuously around the vortex core. 
However, such a kind of vortex-induced S is insufficient, 
at least at low fields, to change the bulk quasi-particle 
DOS N{u) ~ hSnT.J{^ -y6^. + |Afc|2) in that the 
induced S'-component is localized near the vortex core. 
Moreover, the behavior of S does not show any sign of 
increasing sensitivity to the applied field at decreasing 
temperatures. Therefore, we may rule out the role of 
the s-channel in the abnormal thermal conductivity of 
Krishana et ala 

On the other hand, at T < Tsd, T is already broken in 
the uniform phase. Lowest order perturbative treatment 
of S is no longer valid at T <C Tsd ■ S and D are now of 
competing order far from the core with a locked relative 
phase ±7r/2 because of the fact that they are uniform 
there. In other words, symmetry consistency requires 
that S develops a winding component ~ e in addition 
to the e~^^- and e^*^- components, i.e., S ~ ae** -|-&e^*^ -I- 
ce~'^-|-higher order terms near the vortex core. Thus gen- 
erally ISp = a^ + 52 + c2 + 2(a6 + ac) cos 26* + 2&C cos 46*, 
developing a striking two-fold, rather than a four-fold, 
symmetry near the vortex core. The situation is similar 
for the D-component. The relative phase varies again 
near and around the vortex, but is locked to ±7r/2 far 
away from the core. As the two-fold symmetry is gen- 
erated solely from ae*^, the quantity a{r -^ 00) = ztirj 
can also be thought of as the order parameter of the 
two-fold-symmetric vortex structure. Recent numerical 
calculations support the scenario that a structure phase 
transition for the vortex profile convolute with the uni- 



form T-breaking phase transition. t3 Such an exotic vor- 
tex can be possibly observed either by magneto-optical 
spectroscopy, or more possibly by its response to a 'ro- 
tating' applied current J = J{cosu!tx + sinujty). At low 
fields and in a high-quality crystal, there should be a 
resonance in the vortex dissipation at a frequency com- 
parable to the scale of the energy barrier for the vortex 
to rotate by an angle of tt/2 or tt. 

The mixing of the two order parameters S and D gives 
a nontrivial upper critical field, Bc2- In our case, the 
linearized GL equations are: 



2asS + K[2lfS - 



{Ul - nl)D] = ES; 
[Hi - UDS] = ED. 



where we have add an eigen energy term E on the right 
hand side of the equations. E — corresponds to the 
solution of the GL equations. The condition that the 
ground state energy Eg = determines the uppet-critical 
field. This problem has been treated in Refs.M'E3 where 
an implicit solution for Bc2 was obtained. Here we-give 
an explicit solution along the line of Sigrist et al .tj By 
the inspection that [na;,ny] — 2ieB, we can define the 
bosonic operators a and a^ such that [a, a^ = 1, with 



{U, + iUy)/VI^; 



(n, - iiiy)/VI^. 



Here B is the magnetic induction, the fluctuation of 
which at the upper critical field can be neglected. In 
terms of the bosonic operators, the above linearized GL 
equations becomes 

[4:KeB{2n+l) + 2as]S 
+2KeB{aa + a^a^)D = ES, 
2KeB{aa + a)a^)S 
+ [2KeB{2h + 1) + aD]D = ED. 

Here n = a^a is the Landau quantum number. S and 
D can be expanded in terms of the Landau levels \n) as 
{S,D) — X]n(°^n' ^»)l"-)- If there were no coupling be- 
tween S and D, we would obtain the usual result that 
the lowest Landau level n — determines the upper 
critical field. In our case, however, the eigenvalue prob- 
lem involves coupling between all next-nearest Landau 
levels. Expanding the order parameters up to n = 2, 
we can find the ground state energy Eg variationally. 
The upper critical field is obtained by setting Eg = 0, 



6asD + 25]Bo/8 with 



and is Bc2 = [5 - asu + y/^'sD 
Bq = —aD/2eK being the temperature dependent bare 
upper critical field for the pure d.j.2_y2- superconductors. 
In full dimensions, Bq = -~aohc/2eK = ^q/2tt^^ with 
$0 = 2nhc/2e = hc/2e and ^^ = -K/ao- It should 
be pointed out that the above variationaLffcaund state 
is closely similar to that obtained in Refs.EJ't3, where S 
are given by a difference between two gaussian functions. 
This is clear from the fact that S has two nodes in such 
a case, as the Landau state with n = 2 does. Indeed, 



these solutions behave similarly: Bc2 ^ -Bq at T — > To 
(or asD -^ cxd), but it is always larger than i3o, develop- 
ing an upward curvaiuc^in the temperature dependence 
near and below TD.Ilj't3 (It diverges in the unphysical 
region a — > — cxd). This is due to the very effect that the 
sub-dominant channel contributes excess energy lowering 
by adjusting the shape of the vortex. 

Finally, the AbrikosoV|-MQi4ex lattice has also been ad- 
dressed in the literatureJi3il3 The nontrivial four-fold or 
two-fold vortex structures have important impact on the 
lattice structure. Due to the intrinsic anisotropics aris- 
ing fmiH-liS'p they favor generally an oblique vortex lat- 
tice.E3iij Of course, with increasing temperatures, the 
amplitude of S decreases, and one generally expects a 
crossover to the hexagonal vortex lattice near Tq. An- 
other way of seeing this is as follows. Since the vortex 
lattice can be constructed in terms of the (highly degen- 
erated) ground state wave functions for the linearized GL 
equations, it can not be strictly hexagonai-jf the ampli- 
tude of the n — 2 Landau level is finiteJl3 Franz et al 
were also be able to conclude that the orientational an- 
gle of the vortex lattice with respect to the crystal axis, 
a, is no longer arbitrary (as in the case of a single or- 
der parameter), bu|t-|can only take the four possible val- 
ues (X^ ±7r/2,±7rjla as was found in numerical simula- 
tions.Ej In the limit of aso —* oo (or T — > To), however, 
the amplitude for \n = 2) vanishes and we recover the 
usual result for a pure D order parameter. 

Before closing this section, we would like to point out 
that all conclusions drawn for dx2_y-2-\-is superconductors 
are also true for the d^y + is superconductors due to the 
apparent symmetry. 



IV. Dx2_y2 -H/D^ 



SUPERCONDUCTORS 



We consider the case T^ > Td'. Similar to the above 
case, there is also a second order T-breaking phase tran- 
sition. Here the transition temperature T^d' is given by 
a p'D EE a D'/lapl < -1/3 (with an < 0), or Td'd = 
y^T^i/Ti). Below Tddi and in the uniform phase, D' — 
zti-qD with rj = ^— (1 + 3a£)'D)/(3 + a^'D)- However, 
above Tdd' no Z)'-component can be induced by a vortex 
in the absence of a mix-gradient term coupling D' with 
D. Below Tdd', symmetry requires that both compo- 
nents wind in the same manner, e.g., D,D' ^^ e'^. Also 
because of this, the upper critical field is independent of 
ao'D so long as aijiij > —1 (or simply T > 0). In this 
case, D' / D = at the upper critical field, even though 
D' 7^ in the absence of the field. 

In fact, because D and D' are not coupled by gradi- 
ents, the system has two meta-stable states with D' = 
and D' y^ respectively. The stable state is that 
with lower free energy. This is related to the fact that 
the T-breaking phase is frustrated at the vortex core if 
ao'D > ^L While winding simultaneously for D and 
D' increases the kinetic energy of the superfiuid, the 



nonzero D'-component may lower the homogeneous en- 
ergy at low fields. We believe the competing energies 
may drive a transition from the (1^2 _y2 + idxy-'w&Ye vor- 
tex state to a pure dx2_y2-weiYe vortex state at a field 
lower than the upper critical field. By the observation 
that the system is in the dx2-y2 +idxy-Y/a,ve phase at zero 
field but is not at the upper critical field, and that the 
system would be in the d^2_j,2 -|- id^y-wave vortex state 
at all fields when ao'D = — 1, there must be at least one 
field-induced phase transition below Tddi. The transi- 
tion field is zero and Bc2 for T = Tuui and T = 0, re- 
spectively, and should generally increase with decreasing 
temperature. Above the transition line is the region of 
pure dx2^y2 mixed state. Such a transition would predict 
that the quasi-particle excitation gap ( which should be 
proportional to \D'\) decreases with increasing field be- 
low the transition line, and eventually vanishes above the 
transition line, where residual density of states can arise 
from the Donplcr-shift due to the supercurrent around 
the vortex. &0 (The field-induced low-lying quasiparti- 
cle states are possible only for a nodal pairing function) . 
This translates that at T < Tdj^i the thermal conductiv- 
ity should increase rather than decrease with increasing 
field, being consistent with the genesal trend of the data 
at sub-Kelvin temperatures in Ref.cl. But the concrete 
field dependence of the thermal conductivity mayrbe dif- 
ferent. In fact, the mechanism proposed in Ref.cl rules 
out the I?'-component. 



V. STRONG COUPLING EFFECTS: 
PHENOMENOLOGICAL THEORY 

As the electrons in unconventional superconductors are 
strongly correlated because of their d-shell character, we 
now incorporate the strong coupling effects in a pure phe- 
nomenological manner by including the Zeeman energy. 
This energy follows from the interaction between the 
magnetic field and the angular momentum of the Cooper 
pair via strong coupling effects. The underlying mecha- 
nism is simple: Cooper pairs with a definite projection of 
internal orbital momentum will respond to the direction 



of the vortex stmercurrent, which is determine 
magnetic field. t3 The desired energy density isE. 



by the 



-Aireh f d^Rg{R)A^Rx —A 
(X ih{D* D' -c.c), 



R 



(7) 

where h = V x A is the local magnetic induction, g{R) 
is representative of the strong coupling effect, and Ar is 
Cooper-pair wave function in the center-of-mass frame, 
i.e., the inverse Fourier transform of A^ = S+D cos 20^-1- 
D'sin26'fc subject to the weak coupling condition |fc| — 
kp: 

Ar^ t Akf^'^-^dOk 



with Jn{kR) = J cosn0e*'=^™"'^d(/), and 9r being the an- 
gle between R and the a-axis in the copper-oxide plane. 
Here we have adopted the local limit by neglecting the 
spatial variation of the order parameters in evaluating the 
strong coupling term, since Xp ^ ^ with ^ being the co- 
herence length of the superconductor. Note that the final 
result of the strong coupling term does not include the S 
component. This is not accidental but results from sym- 
metry. For convenience we rewrite the non-diamagnetic 
contribution ^) in the form as SkK{2ieh){D*D' — c.c) 
with the only phenomenological dimensionless parameter 
dk. Notice that the Zeeman term can also be rewritten as 
gradient terms 6kK[{UxD)*UyD' ~ (UyDyU^D' + c.c] 
due to the fact that [na;,]!^] = 2ieh. 

To simplify the matter, here we only consider the 
dx2_y2 + idxy superconductor in the presence of strong 
coupling effect. In this case, the new free energy reads. 



F = 



iV(0) 



aD\D\^ + aD'\D'\'' 



n 



^i^i'+^i^r 



\\D\''\D'\^ + \{D*D' + c.c.f 



-K [inZJp + IBD'P + 2ieh6k{D*D' - c.c.)] 



+ (V X Ay/Sn. 
Jn 



(8) 



SJoikR) + {D cos29r + D' sm2dR)J2{kR), 



In the case of au = cto' < and Sk ^ 0, it is known 
that there could be a first-order phase transition from the 
Meissner state to a state with a sudden penetration ofLa. 
finite density of vortices at a lower critical field iJci-t3 
(The transition is of first order since Bd j^ a.t H = 
Hci). However, in the extreme London limit we shall 
adopt (suitable for high-Tc superconductors), we expect 
that both Hci and Bd would be too small to alter our 
subsequent discussion, and will be neglected. 

Let us now discuss the outcome of Eq. (||). To 
be concrete but without loss of generality, we assume 
Tc = Td = lOOX and Td' = lOX. As discussed in 
the previous section, at zero magnetic field, there is a 
second order T-breaking phase transition at a critical 
temperature Tdd' = ^tJ^.ITd « i.X^K. Below Tdd' . 
D' = ±ii]D. Explicitly, \D\'^ = (3 + aD'D)Dl/2 and 
|D'|2 = _(i + 2,aD'D)Dl/2. Here Dl = -aoh- Hence- 
forth we shall normalize the order parameters in units 
of Dq, and the induction field by Bq defined previously. 
They are understood as dimensionless quantities unless 
specified otherwise. Practically for high- Tc superconduc- 
tors, both Dq and Bq saturate at low temperatures (e.g., 
T < mfC): Do - 2Tc - 200K and Bq ~ 5^2(0) ~ AOT- 

iootM 

Following from Eq. (ph, in the presence of the non- 
diamagnetic coupling between the order parameters and 
the magnetic field (5k 7^ 0), a pure dx2_y2 state is sta- 
ble only in zero field above Tddi, while a T-broken 



dx'^-y2zLidxy state is established at all temperatures at fi- 
nite fields, albeit with varying amplitude of D'. We now 
concentrate on the low field regime and in the London 
limit, so that the kinetic energy of the superfluid and the 
inhomogcncity due to vortices can be safely neglected. 
The same methodology has been applied in Rcf.H. At 
T > Tdd' and low fields, the order parameters can be 
treated perturbatively, and are, up to the second order 
in e = SkB (> for ease of presentation) with B as the 
spatial average of the local induction h, 



\D\=2/V3-e\D'\/4, 



(9) 
(10) 



This amounts to a minimum gap on the Fermi surface 

Ao « 2\D'\T, ^ 2€TdTd'd/{T - Td'd) 
(xSkB/{T-TD'D) (11) 

at e ^ 1 and T > Td'd (not very close to Td'd)- Ne- 
glecting for the moment the suppression due to quasi- 
particle scattering at the vortices, we observe that the 
thermal conductivity Kg at T > Tb'd in the presence of 
a minimum gap Aq is given by. 



Ke{T,B)/K,{T,Q) ^ (l/T)exp(-Ao/T) 
- (1/T) exp[-const. x 5kB/{T ~ Td'd)T\. 



(12) 



It should be appreciated that this exponential decay in 
B is increasingly sharpened by T[T — Td'd) instead of 
only T. In Eq. ^ a power law Ke(T,0) ex T^ jg ^ged, 
which follows from the fact that above Td'd the zero- field 
state is a pure d.j.2_y2 state. Further we have implicitly 
neglected the change in \D\, which is merely of the second 
order in e from Eq. (g). The exponential decay develops 
a crossover at Aq ^ T, or. 



Bk<xTiT~TD'D) 



(13) 



which indeed resenibles a scaling law Bk oc T^ addressed 
by Krishana et ala As an example, we plot \D'\ and 
Ke{T,B)/Ke{T,0) in this context in Figs. 1(a) and (b), 
respectively. Here \D'\ is calculated exactly from the 
optimization of the free energy subject to the neglect- 
ing of kinetic energy and field energy, which should be 
reliable at low fields. Defining a criterion for the kink 
transition (more precisely, cross-over), the parabolic- 
like contour lines unambiguously reproduce the scaling 
law for the kink field.Q The sharpening feature of the 
crossover is also clear. Including the additional sup- 
pression of Ke from quasi-particle scattering by vortices 
would render even better qualitative agreement with 
the experimental results, but we shall not go into fur- 
ther details. Interestingly, Eq. ( [l3| ) predicts a down- 
ward curvature of the Bk(T) curve in a log-log plot 
on the low temperature, side, which seems to be the 
case in Fig. 4(b) of RefJj. We have also performed ex- 
act calculation of the full set of GL equations ensuing 



from Eq.(||) to find that the results in the low field 
regime remains essentially unchanged, while at inter- 
mediate fields, \D'\ saturates a while and then dimin- 
ishes gradually together with \D\ near the upper critical 
field Bc2- The latter is given by, in dimensionless form, 

Bc2 = 1 - ao'D + ^/il + ao'Df-^aD'oSk^ /2(1 - 
(5fc2). At Bc2 the ratio \D'\/\D\ is given by \D'\/\D\ ^ 
5kBc2/{oLD'D + Bc2)- This means that T-breaking is re- 
tained up to the upper critical field. However, it is clear, 
in conjunction with Fig. 1(a), that the strong coupling 
effect vanishes near the critical temperature T^ (where 
ctD'D -^ oo), a necessary ingredient for the theory to be 
compatible with the general properties of high-Tj. super- 
conductors at temperatures near Tc. Also interestingly, 
in the context of this theory, (i) the power law Kg oc T^ is 
invalid at finite fields above To'd, which seeins to be con- 
sistent with the experimental data in Refill, and would 
be violated below To'd even in the absence of magnetic 
field; (ii) Kg increases with increasing magnetic field be- 
low To'D- On the other hand, the theory does not involve 
a phase transition for the development of a plateau in Ke, 
because there is no further symmetry breaking at finite 
fields in our theory. 








T(K) 



FIG. 1. Gray-scale contour plots of (a) \D'\ and (b) 
Ke(r,B)/«:e(T,0) fa (l/T) exp(-Ao/r) as functions of B and 
T. Here Sk = 1/2 and Aq = 2\D'\Tc with Tc = WOK. 
Tod' ^ 3.16A" is the zero-field T-breaking transition tem- 
perature. See the text for details. 



VI. SUMMARY 

We have developed a GL theory with respect to the 
three relevant singlet pairing channels in unconventional 
superconductors. The case of d -I- is superconductors 
has been discussed. The novel two-fold symmetric struc- 
ture of vortices below the T-breaking temperature is pre- 
dicted. The case of d^^-y^ -\- id^y superconductors has 
been addressed in detail. In the absence of strong cou- 



pling effects, we find tfiat tfie sub-dominant d^y-wave or- 
der parameter can not be induced by vortices above the 
T-breaking temperature, in contrast to tfie case in a d+is 
superconductor. Below the T-breaking temperature, we 
predict a novel field-induced first-order phase transition 
from the d^^-y^ + id^y vortex state to the dj.2_y2 vortex 
state at higher fields. By including the strong coupling 
effect phenomenologically, the d^2_y2 +idxy state can be 
induced by the magnetic field even above the zero-field 
T-breaking temperature. Near and above this transition 
temperature, the induced da;y-wave order parameter may 
be compared with some results ior the abnormal ther- 
mal conductivity reported in RefMI. Although our theory 
could give the general trend that the thermal conduc- 
tivity would increase with increasing field below the T- 
breaking temperature, it can not result in a pure d^i^yi 
state near the zero temperature, which was claimed to be 
essential-to explain the sub-Kelvin thermal conductivity 
in Ref.d'EZI. If that is the case, we argue that there are two 
possibilities for the material used in RefJj: either some 
significant change occurs for the d^jy-channel near zero 
temperature such that the da;j,-wave order parameter is 
suppressed drastically or no (ia;j,-channel is present at all. 

It should be pointed out that the s- and dxy- channels 
are only two, but not all of possible sub-dominant singlet 
pairing channels in a d^i^-ys -dominant superconductor. 
We have neglected, e.g., the kxky{k^ — fc^)-channcl (in 
the A2g representation), which would appear in the case 
of a more general pairing interaction, or in a higher-order 
expansion of the gap equation. lj This order parameter 
could also be induced by spatial inhomogenieties of the 
dx2-y2 order parameter, and would indirectly induce a 
dxy order parameter. However, odd-parity order param- 
eters can not be induced in singlet-pairing superconduc- 
tors. 

In this paper, we have only performed accessible an- 
alytical as well as qualitative discussions on the theory. 
More exact and numerical results are awaited. The the- 
ory is also highly useful for studying-Lhe vortex dynamics 
of unconventional superconductors. t3 
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